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Abstract: We analyze the contributions to rare kaon decays mediated by fla¬ 
vor-changing Z-penguin diagrams in a generic low-energy supersymmetric extension 
of the Standard Model. In order to perform a model-independent analysis we expand 
the squark mass matrices around the diagonal, following the so called mass-insertion 
approximation. We argue that in the present case it is necessary to go up to the 
second order in this expansion to take into account all possible large effects. The 
current bounds on such second-order term, which was neglected in previous analyses, 
are discussed in detail and the corresponding upper bounds for the rare kaon decay 
rates are derived. As a result, we show that supersymmetric effects could lead to 
large enhancements of K ttuu and Kl 7T^e'^e~ branching ratios. 
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1. Introduction 

Flavor-changing neutral-current (FCNC) processes provide a powerful tool for indi¬ 
rect searches of New Physics. This is particularly true in the framework of low energy 
supersymmetry |]I[, which represents one of the most interesting extensions of the 
Standard Model (SM). The large number of new particles carrying flavor quantum 
numbers, present in this context, would naturally lead to sizable effects in FCNC 
transitions [Q, |]. 

At the one-loop level, supersymmetric contributions to FCNC amplitudes can 
be classihed into three groups, according to the virtual particles inside the loop: 
i) Higgs/W-quarks, ii) gluino-squarks and hi) chargino/neutralino-squarks. The 
hrst group contains the SM contributions as a particular subgroup, whereas ii) and 
hi) represent genuine supersymmetric effects. Among them, gluino-squark transi¬ 
tions have been widely discussed in the literature i, I, I, i and are expected to 
produce the dominant non-SM effect in AF = 2 processes. This is conhrmed by 
the analysis of mixing with the inclusion of gluino-squark contributions, 

which provides severe constraints on supersymmetric models [^, p. The effect 
of chargino/neutralino-squark diagrams is usually neglected in the analysis of such 
processes. 
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A different situation occurs in AF = 1 transitions mediated by Z-penguin di¬ 
agrams, which are particularly relevant to rare kaon decays, like K ttvV. As 
recently discussed in ii, the dominant supersymmetric contribution to these pro¬ 
cesses is given by chargino-up-squarks diagrams. This is because the Zqiqj effective 
vertex is necessarily proportional to S'17(2)2,-breaking couplings that, in supersym¬ 
metric models, are provided by qi — qn, qi — Qr and wino-higgsino mixing. Since 
the qi — qR mixing in the down sector is suppressed by the small down-type Yukawa 
couplings, the effect of gluino and neutralino diagrams is necessarily small. On the 
other hand, the large Yukawa coupling of the top leads to potentially large effects 
in diagrams involving up-type quarks or squarks. Indeed the (d, s)l — mixing, 
already present in the SM, is responsible for the mf enhancement of the Higgs/W- 
quark contribution to the Zsd effective vertex. Analogously, it is natural to expect 
a large effect due to the (d, s)l — tn mixing in diagrams involving charginos and 
up-type squarks. 

This effect has been already noted by Buras, Romanino and Silvestrini in the 
calculation of the supersymmetric contributions to K ^ ttz/z/ [Q. However, this 
calculation has been performed in the single mass-insertion approximation, where 
only terms with at most one off-diagonal element of the squark mass matrix are 
considered. We believe that this approximation is not sufficient to fully account for 
possible large effects in the present case. Indeed, in order to provide the necessary 
SU{2)l breaking, at least two mass-mixing terms are necessary, either from the 
squark sector {ul—ur) or from the chargino sector (wino-higgsino). Since both these 
couplings vanish as 0{mw/Ms) in the limit of a heavy supersymmetry-breaking 
scale {Ms), we consider more appropriate to expand in both of them up to the 
second order. One could argue that wino-higgsino mixing is not suppressed by 
the off-diagonal flavor structure. However, the hierarchy of the Yukawa couplings 
implies that terms with a single wino-higgsino mixing always appear together with 
suppressed CKM factors. As a result, it is reasonable to expect that terms with a 
double LR mass insertion and without any CKM suppression are at least of the same 
order as those generated by a single LR mass insertion together with wino-higgsino 
mixing. 

In the present paper we present a complete discussion of the supersymmetric 
contributions to the Zsd amplitude, beyond the single mass-insertion approxima¬ 
tion. We hnd that the contribution generated by a double LR mass insertion in 
the up-squark sector, which was neglected in previous analyses, yields a potentially 
large effect. Employing the notation of [§], we can naively say that in this case the 
factor of the SM amplitude gets replaced by In¬ 
terestingly, this kind of mechanism is only weakly constrained by mixing 

and can provide a sizable enhancement (up to two orders of magnitude) to rare decay 
widths. 

The paper is organized as follows. In section 2 we discuss the supersymmetric 
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contributions to the Zsd amplitude, with particular attention to the hierarchy of 
the various terms. The role of box diagrams in AF = 1 transitions is also briefly 
analyzed. In section 3 we discuss theoretical and phenomenological bounds on the 
up-type LR couplings. In section 4 we analyze the possible enhancements of rare kaon 
decays rates driven by these supersymmetric effects. The results are summarized in 
the conclusions. 

2. The Zsd effective vertex 

The amplitude we are interested in here is the one-loop FC effective coupling of 
the Z boson to down-type quarks, in the limit of vanishing external masses and 
momenta. As already emphasized in [^ |^, the SU{2)l x U{1)y —> U{l)e.m. gauge 
structure implies that this coupling proceeds through symmetry-breaking terms and 
involves only left-handed quarks. Thus it can be generally described by introducing 
the effective Lagrangian 



( 2 . 1 ) 


where Wds is a complex dimensionless coupling. 

In our conventions, the SM contribution of top-quark penguin diagrams, evalu¬ 
ated in the’t Hooft-Feynman gauge, leads to 


ir/y = \c(x ,„), 


( 2 . 2 ) 


where A* = V^gVtdi Vij are the CKM matrix elements 0 and Xtw = The 

loop function C{x), originally computed in [^, can be found in the appendix. We 
recall that C{x) —>■ x/S for large x. 

In the minimal supersymmetric extension of the SM, which requires two Higgs 
doublets, the contribution of penguin diagrams with the exchange of charged Higgs 
and top-quark is aligned with the SM one (i.e. is proportional to A*). Denoting as 
usual by tan /3 the ratio of the two Higgs vacuum expectation values [Q , we find 



where now XtH = Similarly to the SM case, also H{x) —>• x/S for large 


X (the full expression of H{x) is given in the appendix). The sum of (|2.2|) and 
(|2.3|) complete the hrst class of contributions outlined in the introduction, namely 
the Higgs/W-quark diagrams. To analyze the genuine supersymmetric effects, and 
particularly those generated by chargino-squark exchange, we hrst need to discuss 
shortly the structure of the supersymmetric mass matrices. 
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In the basis of the electroweak eigenstates, wino and higgsino, the chargino mass 
matrix is given by 


M^ = 


M 2 \/2mw sin f3 
\/2mw cos (3 /i 


(2.4) 


where the index 1 of both rows and colnmns refers to the wino state. Following the 
standard notation [^j, here fi denotes the Higgs quadratic coupling and M 2 the soft 
supersymmetry-breaking wino mass. To dehne the mass eigenstates we introduce 
the unitary matrices U and V which diagonalize My. 




(2.5) 


As can be noticed, the off-diagonal entries of My are 0{mw), whereas M 2 is 0{Ms). 
In the limit where m\Y/M 2 is a small parameter we can perform a perturbative 
diagonalization of My around its diagonal elements, or, correspondingly, an expansion 
of U and V around the identity matrix. 

In the squark sector we have 6x6 matrices which mix the three families of 
left-handed and right-handed squarks. A convenient basis for our calculation is the 
basis where the d\ — v3i^ — Xn coupling is flavor diagonal and the — 12]^ — Xn one is 
ruled by the CKM matrix (see [§ for a more detailed description). In this case, the 
up-squark mass matrix is given by the Hermitian matrix 


M 


2 _ 
u — 


( {^u)dLdL {^u)dLUR \ 
\{Ml)uRdR {MI)urUr) 


( 2 . 6 ) 


where the subscript di (which runs over three values) indicates the combination of 
left-handed up-type squarks which appear in the diagonal couplings — Xn- 

On the other hand, the index ur denotes the combination of right-handed up-type 
squarks which appear in the d\ — vdp(^ — Xn vertices ruled by the CKM matrix-element 
Vij. Since is Hermitian we need to introduce only one unitary matrix, id, to 
diagonalize it 

HMlH^ = . (2.7) 

Also for M-ij the off-diagonal elements are expected to be small and the perturbative 
diagonalization is well justihed 0. 

We are now ready to evaluate the contribution of the chargino-squark penguin 
diagrams in Fig. |^. The full result before any mass expansion is quite simple and is 
given by^ 

= ^A'‘,A‘,F,U, , ( 2 . 8 ) 

^ The sum over the repeated indices i and j (running from 1 to 2), Z and k (running from 1 to 
6), and (running over the three values d^, sl and ) is understood. 
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where 


4 = Hu, V}j - giVuHu, V^, , (2.9) 

At, = HlJ',, - gKHltVa . ( 2 . 10 ) 

^jilk ^Ik ^iS^iki^jk) 5i}^ yj^ik^jkj i^^iki ^jk) 

^ij /c(Xj^, (C;fc) . (2.11) 

For simplicity the effect of all the Yukawa couplings but gt = rrit/{\/2mw sin (3) has 
been neglected. Analogous to the previous cases, the variables Xij denote ratios of 
squared masses (e.g. Xik = M‘^JM‘1^) and the functions k{x,y) and j{x,y) [§ can 
be found in the appendix. 



Figure 1: Chargino-up-squark penguin diagrams contributing to the Zsd effective vertex 
(diagrams involving self-energy corrections to the external legs are not explicitly shown). 

The product of Aji and A|^ in ( |2.8|) generate four independent terms, proportional 
to g^Xt, gtVtd, 9tVt*s 1, respectively, which correspond to the so-called RR, LR, 
RL and LL contributions in the notation of 0,1- We proceed to analyze these terms 
separately. 


1. The RR contribution is generated by a Yukawa-type interaction in both quark- 
squark-chargino vertices of Fig. |^. This term is the only one which survives in 
the limit of diagonal Aifj, i.e. to the lowest order in the perturbative expansion 
of H around unity. In this limit is given by 


^L\rr — g9t^tV2j 


Uij ^/Xit^Xjtnj {xit^ , Xjtn ) 


R: 


i2 


( 2 . 12 ) 


As anticipated in the introduction, the SU{2)l breaking of the Zsd vertex 
requires at least two mass-mixing terms, either from the squark sector or from 
the chargino sector. In (|2.12|) the absence of the former mechanism implies 
a double wino-higgsino mixing, as can be easily checked by the mismatch of 
V and f/ indices. Thus W^J^rr is parametrically suppressed by 0{rn^/M 2 ) 
and aligned in phase with respect to . To get a feeling of the numerical 
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factors, note that /c(l, 1) = 3/2 and j(l, 1) = 1/2. We then conclnde that this 
contribution cannot provide a sizable effect, particularly in the limit of a heavy 
supersymmetry-breaking scale. 

Considering higher orders in the perturbative expansion of H, one can easily 
check that there is no contribution to at the first order. At the second 

order it is possible to generate a non-vanishing contribution and also to avoid 
the wino-higgsino mixing. However, the unavoidable factor Xt makes W^g\iiR 
always not particularly interesting with respect to . 


2. The LR and RL terms are originated by a Yukawa-type interaction in one of 
the two quark-squark-chargino vertices of Fig. and a gauge-type interaction 
in the other. As can be easily understood, this implies that W^^Ilr and 
are at least of hrst order in both wino-higgsino and qi — Qr mixing. Performing 
explicitly the expansion of H up to the hrst order, as discussed in the appendix, 
we hnd 


= -latVu 


M? 

UL 


5ijh{Xiu]^1 1) ‘^Ui\Uyj (X iuLi Xtuui)\Vii , (2.13) 


where Mu^ indicates the average mass of the approximate left-handed up 
squarks. The term can be obtained from (p.l3|) with the substitu¬ 
tion Vtd{Ml)s^tn Vt*siMu)tRdL and with the exchange 1 2 in the indices 

of the V matrices outside the square brackets. 


The presence of a single CKM matrix element in (p.l3|) leads to potentially 
large effects: the missing factor V)* is replaced by where 


(ilnU = , (2.14) 

and the ratio {5^R)tslVts can be larger than one |^. However, this enhancement 
is partially compensated by the 0{mw/M 2 ) suppression induced by the wino- 
higgsino mixing and the total effect is not very large. Indeed the phenomeno¬ 
logical bounds on dictated mainly by 6 —>• sy, become weaker for large 

supersymmetric masses, when the wino-higgsino suppression gets stronger. As 
a result, can be at most as large as [||. Similar comments apply 

also to (see the next section for a more detailed discussion about limits 

on i6^R)ts and 

It is interesting to note how the LR and RL terms, which arise only at hrst 
order in the expansion of H, are potentially larger and not aligned in phase 
with respect to the lowest order contribution, provided by This is 

clearly a consequence of the disappearance of one of the two CKM factors. For 
this reason, it is natural to expect that terms arising at the second order in 
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the mass-insertion expansion and without any CKM suppression could be even 
bigger. 


3. The LL term is originated by a double gauge-type interaction in the quark- 
squark-chargino vertices of Fig. |^. Similarly to the RR case, this implies a 
second-order mixing either in the chargino sector or in the — Qr sector. 
However, contrary to the RR case, there is no contribution to the leading- 
order in the expansion of H. The hrst term in this expansion arises to the hrst 
order and is given by 




1 i^u)sLdL 


M? 

ul 




ij 


V,2V^^k{x.. 




, 1 ) 


-2Ui2U2j {x 


; Xjuj ^) 


,l)Fi. (2.15) 


As can be noted, this term involves a double wino-higgsino mixing, which 
provides the necessary SU{2)l breaking, and a hrst-order mixing among left- 
handed squarks. Thus, even if apparently enhanced by the absence of any 
CKM factor, is strongly suppressed in the limit of heavy supersymme¬ 

try-breaking scale. Moreover, the SU{2)i invariance of the soft-breaking terms 
relates {M^)s^dL fo {My)s^dL 0) which is strongly constrained by 
mixing. As a result, turns out to be always smaller than p. 

A different scenario occurs if we consider the contribution to 1 T^|ll which 
survives in the absence of wino-higgsino mixing. In this case one has to go at 
least to the second order in the expansion of FT, and only terms with a double 
LR mixing survive. The lowest-order result in this limit is simply given by 


^/X^2 ^ 1 i^u)sLqRi^u)qRdL 

ds I LL '* 


m 




XI 


Q i^LR}qs i^LR}qd ^uj^l Xq^i) , 


(2.16) 


where the function l{x,y,z), normalized to /(1,1,1) = —1/12, is reported in 
the appendix. Contrary to the cases of WdsliR ^^nd discussed 

previously, there is no explicit suppression in fhe limit of heavy 

superpartners. Actually the upper bounds on the factor go to 

zero in this limit, as we will see in the next section. However, for q = t there 
is room enough to produce sizable effects (in close analogy to the Xt factor in 
the SM case) even for Mr ~ 1 TeV. In this case could be substantially 

larger than providing sizable enhancements to rare kaon decay rates. 


Using the effective Lagrangian (|2.1|) we can easily calculate the effects of the Zsd 
penguins discussed above in various processes. In the case of K ^ nuu decays we 
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( 2 . 17 ) 


find that the contribution generated by (^?T|) to the X function, dehned by 
G p a 


^e// — 


v/2 27rsin^0 


w 


XtX S 7 ^(l - 75)d i/z 7 ^(l - 75 )z/z + h.c., 


is given by 

^Zsd = Wds/Xt ■ (2.18) 

Comparing our results in (|2.3|) , ( p.l2|) , (|2.13|) and (|2.15|) with those reported in the 
appendix of we hnd a perfect agreement but for an overall factor —1/2 which is 
a mere misprint.^ 

In principle, in the case of iC — ttz/P decays also supersymmetric box diagrams 
could provide sizable effects, as it happens for instance in the SM case [jTT[ . However, 
the contribution of chargino-up-squark box diagrams to X turns out to be always 
suppressed by a factor besides possible wino-higgsino mixing.^ In a generic 

expansion in powers of off-diagonal mass terms, denoted by e, the box contribution 
to X starts at 0{e^), whereas the penguin one at 0{e^). Thus, in general we agree 
with the statement of Nir and Worah 0] that penguin contributions provide the 
dominant effect. Only in the case of the terms proportional to when the 

penguin contribution is suppressed and starts at C>(e^), the corresponding box term 
turns out to be competing ||^. However, as long as we are interested only in possible 
large effects this is not a relevant case. 

Similar arguments apply also to other processes where the effective Zsd vertex 
can contribute, like K —> and K —> decays. Hence, in a minimal 

supersymmetric extension of the SM with generic flavor sector, and particularly in the 
limit of a heavy supersymmetry-breaking scale, we consider it a good approximation 
to encode the dominant non-SM effects to these processes via the Lagrangian m- 
A similar approach was considered by Nir and Silverman in a different context [|T2 
the coupling Wds in (|2.1|) is related to the Uds of Nir and Silverman by 


Ud.= 


a 


TTSin 


20 


-IT, 


ds 


W 


(2.19) 


3. Bounds on the couplings. 

In the previous section we have argued that the (LR)^ term that appears at second 
order in the mass-insertion expansion, may give the largest enhancement to the 
Zsd effective vertex with respect to the SM contribution. In the present section 
we will analyze in detail the bounds we can put on this term, considering both 
phenomenological information, and purely theoretical constraints. 

^ This misprint does not affect the numerical results of |^. We thank L. Silvestrini for clarifying 
this point. 

^ The contribution of the charged-Higgs box diagrams is clearly negligible because of the small 
lepton Yukawa couplings. 

















3.1 Vacuum—stability bounds. 

Before analyzing the bonnds coming from phenomenology, we discuss an interesting 
result obtained by Casas and Dimopoulos [jT^, who have shown that bounds on the 
off-diagonal LR entries of the squark mass matrices can be derived also from the 
requirement that the standard vacuum of the theory be stable. In particular they 
require that there are no charge and color breaking minima (CCB bounds), nor di¬ 
rections in which the potential is unbounded from below (UFB bounds). Obviously, 
these bounds have to be satished by any model (and interestingly enough, are gener¬ 
ally not satished). The only way to avoid or at least to soften these constraints is to 
assume that we live in a sufficiently long-lived metastable vacuum. However, to be 
more conservative, we will not take into account this possibility. The consequence of 
the stability bounds for the matrix elements of our interest can be stated in a very 
simple manner: 




LR)ij 


< m 


2M| + Mf 




Mi 


{k = max(q j)) 


(3.1) 


where M| and M~ denotes the average masses of up-squarks and sleptons, whereas 
niu^ indicates the mass of the Uk quark. The factor provides a very stringent 
suppression unless one of the two generation indices {i and j) is equal to 3. For this 
reason it is a good approximation to replace the sum Y.q{.^^R)*qs {^LR)qd in ( P-16|) with 
the product {^LR)td- In this case the bound (|3.1|) can be roughly expressed 

in the following form 




LR)tsi^. 


LR)td 


< 


3mf 

W 


(3.2) 


where with Ms we have indicated a typical supersymmetric scale. Actually the 
bound ( p.l|) corresponds to the UFB constraint, but as long as we consider almost 
degenerate supersymmetric particles CCB and UFB bounds are essentially equivalent 

0- 

At this point it is useful to make a hrst estimate of the possible enhancement 
induced by the (LR)^ mass insertion in the Zsd vertex. Comparing (|2.2|) and (|2.16|) , 
in the limit x^w 3> 1 and assuming almost degenerate supersymmetric particles, 
leads to 




ivf," 



i^LR)ts i^LR)td 


12 A 

12 xtw At V Ms J 

where the last inequality has been obtained imposing the bound ( |3.2|) . As can be 
noticed, though stringent the model-independent constraint leaves enough room for 
a large enhancement, even for Ms as large as 1 TeV. 


3.2 Box AR = 2. 

A term with two LR mass insertions appears in the box diagram (containing charginos 
and squarks) contributing to K^—K^ mixing. In this case, however, this term appears 
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only at a subleading level. The complete expression for the contribution of the box 
diagram to the effective Hamiltonian for AS* = 2 is 


^AS=2 = 0 ^ k{Xik,Xjk,Xik){sLl''dL){sLlM • (3.4) 

We may now expand the mass matrices around their diagonal part to the desired 
order. Considering only terms without chargino mixing, we get 


-^ik^jk^jl^il k{Xik-i Xjk, Xlk) — 


1 

10 


m,ur - + 


u 
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{i^LR)ts{^LR)td + ...)+ 0{X‘f) . 


(3.5) 


To obtain this result we have not only applied the formulae for the perturbative 
diagonalization of the mass matrices (that we give in appendix), but have also taken 
the limit where all superpartners have approximately the same mass {xki = 1 for all 
fc’s and i’s). If we now use the experimental information on Auik = 3.5 x 10“^^, and 
require that the contribution of the term with two LR insertions in (^)^ does not 
exceed the experimental value, we get 



( Ms \ 

VhOOGeVy 


(3.6) 


We remark that this limit is derived using the quadratic term in (|3.5|) , as the linear 
one is multiplied by Xt which suppresses its contribution strongly. Similarly, we 
have not considered the bounds that could be obtained on the single and 

i^LR)td couplings, which always appear suppressed both by CKM factors and chargino 
mixing. Of course this limit is rather generous, as one would expect the first two 
terms in ( p.5|) to be responsible for the main part of the effect. On the other hand, 
until we will be able to get some independent information on the hrst two terms in 
the expansion (and on their signs too) this is the best we can get from this quantity. 

If we now look at the imaginary part of the same matrix element, and consider 
the experimental information on Re(e), we can get a bound on the imaginary part 
of the (LR)^ term squared: 



/ Ms \ 

UoOGeVy 


(3.7) 


3.3 Limits from B and D physics. 

Buras, Romanino and Silvestrini have analyzed the bounds on various mass in¬ 
sertions coming from R-meson phenomenology. From the chargino contribution to 

^ Evaluating this with the approximation {K^\{sL^^dL){sL'^f_idL)\K°) = 1/3 rriK Pk- 
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Bd — Bd mixing they get 


(5 


U ' 
LR, 


dt 


< 0.1 X 


f Ms \ 

UoOGeV/ 


(3.8) 


A bound on the other matrix element of our interest was earlier obtained by Misiak, 
Pokorski and Rosiek analyzing the chargino contribution to 6 —sy [^: 




Li?) St 


< 3 X 


/ Ms 
VhOOGeVy 


(3.9) 


In principle a limit on Be[{{6^j^)lg{Sspi)tdY], similar to the one in (^.6[), could 
be obtained from the analysis of the gluino-up-squark box diagram contributing to 
D^ — D^ mixing. Note, however, that this bound is very different from those discussed 
above since it can be made arbitrarily small in the limit of a heavy gluino mass. 
Assuming gluino and chargino approximately degenerate, the constraint obtained by 
— Z)° mixing is essentially equivalent to the one. Indeed the ((y'strong/5')^ 

enhancement of the gluino box diagram with respect to the chargino one is almost 
completely compensated by the less stringent experimental constraint on Am^i with 
respect to ArriK-^ 


3.4 Bounds on the Zsd vertex. 


As anticipated in the previous section, the Zsd effective vertex contributes to various 
rare kaon transitions. Some of them have been observed, whereas stringent experi¬ 
mental limits exist on the others: we can therefore use these informations to derive 
bounds on the (LR)^ term which we are now analyzing. These bounds are best 
expressed in terms of the coupling Wds introduced in ( |2.1| ). A similar analysis has 
been already made by Grossman and Nir |]^, using exactly the same language of an 
effective Zsd coupling (but using the Uds of 0)- Following and partially updating 
(and correcting) their results, we find 


1. from the process Ki 




15 : 


|Re(IRds)| <2.2 X 10 


-3 . 


(3.10) 


2. from® B{K~^ —> < 2 x 10 ® [|I^ : 

|1V*| < 3.6 X 10-=> ; 


(3.11) 


3. from the measurement of e: 

\Re{Wds)lm{Wds)\ < 1.1 x 10"® . (3.12) 

® We are grateful to M. Worah for a clarifying discussion about this point. 

® Note that the corresponding bound in |Q was larger due to missing factor six in their Eq. (13). 
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In principle similar bonnds could be obtained from e'/e and B{Kl 7r°e+e“). 
However, in both cases the large theoretical uncertainties and the poor experimental 
information lead to weaker constraints. To translate the results (|3.10| - |3T^) into 
bounds for the (LR)^ term of our interest we use the relation 

^dslh = '^i^LR)tsi^LR)td , 


derived from ( p.l6|) in the limit of degenerate supersymmetric particles. We then 
obtain 


\iS^LR)WUtd 


< 0.21 , 

< 0.35 , 

< 0.1 . 


(3.13) 


3.5 Summary of the bounds. 

The two limits derived from the analysis of the AS = 2 box diagram can be written 
as follows: 


1 2 


) L () irf ) I™ ( () L () 
those obtained from B physics lead to 


i^LR)tsi.^LR)td 


< 0.3 X 


( 


whereas the model-independent one is given by 


i^LR)tsi^LR)td 


< 0.3 X 


<2.6 X 10“2 X 

/ Ms 

UoOGeV 

<1.1 X 10“^ X 

/ Ms 

UoOGeV 

CO 


500GeVy ’ 


500GeVy 


Ms ) • 



) .(3.14) 
) ;(3.15) 

(3.16) 

(3.17) 


Finally, the ‘scale-independent’ limits derived from the phenomenological analysis 
of the Zsd vertex are 


Re ((3" 


LRJts i^LR)td 


i^LR)tsi^LR)td 


U 


< 0.21 , 
< 0.35 , 


(3.18) 


where we have skipped the bound on the product of real and imaginary part, which 
is clearly negligible with respect to the one in (|3.15|) . 

A summary of the various bounds is displayed in Fig. ^ for the sample value 
Ms = 500 GeV. From the hgure it is clear that the bound in (|3.15|) is by far the 
most stringent one. This implies that, if we assume that real and imaginary parts of 
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Figure 2: Summary of the bounds on the (LR)^ coupling At, as defined in Eq. (3.19) 


i^LR)tsi^LR)td are of the same order, these are 0(10“^). On the contrary, if one of the 
two is zero the other can be 0(10“^). The maximum value allowed for the real part 
setting the imaginary part to zero, or vice versa, is 0.16 as dictated by the bound 
(|3.14|) . Playing around with the Ms dependence of these bounds one can hnd that 
the maximum value allowed for either the imaginary or real part (when the other is 
set to zero) can grow up to 0.2 for Ms = 600 GeV, again a bound dictated by ( |3.14|) . 
Above this value of Ms the model-independent limit on the modulus becomes more 
stringent. Notice in fact that the model-independent limit on {Sis.)ts{^LR)td and 
the one coming from S-physics have opposite dependence on the average mass of the 
superpartners. So that for Ms < 500 GeV it is the S-physics one which dominates, 
whereas above 500 GeV the model-independent one takes over. 

In conclusion, for Ms ^ 600 GeV we can just consider the two bounds in ( p.l5|) 
and (p.l7|) , as all the others will be automatically satished. If we dehne 


the bound we have to satisfy 
|At| < min 


\ — {^LR)tsi^LR)td , 


for Ms ^ 600 GeV is 

2 X 10-2 


/ 600GeV \^ 2 x IQ-^ / Ms A 

' y Ms ) ’ UoOGevi 


(3.19) 


(3.20) 


whereas the phase 9t is unbounded. 


13 






























At this point one conld argne whether a reasonable low-energy snpersymmetric 
model conld satnrate this bonnd. It is beyond the scope of this paper to analyze 
any model in detail. However, we recall that in generic snperstring scenarios the 
so-called A terms, responsible for the LR entries of the mass matrices, are expected 
to be 0{Ms) []^. This wonld imply ~ 0{mw ■ Ms). Thus in general it is 

not unnatural to consider models where the bound ( |3.20|) is saturated (see e.g. the 
discussion at the end of Ref. 
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4. Phenomenological consequences of a large Zsd effective cou¬ 
pling 


The most clear signature of an enhancement in the Zsd effective vertex could be 
found in K ^ vrz/P decays. Within the SM these transitions can be described by 
means of the Hamiltonian ( p.l7|) , with the X function given by 

XsM = Xtixtw) + ^Pc, (4.1) 

where Xt{xtw) — 1-5 is generated by the dominant top-quark contribution (summing 
penguin and box diagrams) and Pc = 0.40±0.06 is due to the charm loops (as usual A 
denotes the Cabibbo angle and Ac = V*gVcd) thus |A^Ac/Af| ~ 0{1)). The branching 
ratios of and Kl modes can be expressed in terms of the X function as 



BR{K+ TT+iyp) 

BR{Kl vrVz/) 


K+ 


\tX 

X 


(. \tX 


(4.2) 

(4.3) 


where = 4.11 x 10“^^ and kl = 1.80 x 10“^° [|^. For a numerical estimate we 
recall that A = 0.22, \Xt\ ~ 3 x 10“^ and ImAi ~ |At|/3 [^. 

In extensions of the SM where the main new-physics effects can be encoded via 
the effective coupling Wds, we should add to Xsm the Xzm function defined in (|2.18|) . 
Thus if we add to the SM contribution the dominant supersymmetric effect, provided 
by the (LR)^ terms, we find 


1A A^A 

Xtot ^qnl) P Xf^ (^X^\y) H : Pr 


8 At 

1 At , A^Ac 

qFV (^iw) H-7— 

96 At At 


At 


Pc 


(4.4) 


where At has been defined in (|3.19|) and the second line of (^4.4|) is obtained in the 
limit of almost degenerate superpartners. Given the constraints on |At| reported in 
(|3.20|) it is clear that large enhancements with respect to the SM case are possible. 
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In the rate of the charged mode one can gain up to an order of magnitude if Mg 
is around 600 GeV, where the effect is maximum. In the Kg case a crucial role is 
played by the new phase 9t: if 9t ~ 90° a huge enhancement (up to two order of 
magnitudes) is possible for a wide range of Mg. In Table we summarize the upper 
bounds on the two modes for Mg ~ 0.6 TeV and Mg ~ 1 TeV. 


decay mode 

maximum SUSY branching ratio 
Mg ~ 0.6 TeV Mg ~ 1 TeV 

SM branching ratio 


1 X 10“^ 

4 X 10-1° 

(9.1 ±3.8) X 10-11 19 

Kg —>■ 

4 X 10“^ 

6 X 10-1° 

(2.8 ± 1.7) X 10-11 19 

Kg —>• 7r°e+e“ 

6 X 10-1° 

1 X 10-1° 

< few X 10-11 19, 20 


Table 1: Approximate upper bounds for the branching ratios of A' —> ttuP and 

^ 7r‘^e"*'e“ decays within the low-energy supersymmetric scenario discussed in the 
text, compared to the SM expectations. 


Related modes which could allow to detect an enhancement in the Zsd effective 
amplitude are the K —>■ decays. In the charged channel the long-distance 

process —>■ 7r’''7* —>■ is by far dominant, hiding the contribution of the 

Zsd transition. However, the single-photon exchange amplitude is forbidden hy CP 
invariance in the Kg —>■ mode, which is therefore more sensitive to short- 

for a recent discussion about these decays), 
and Kl ^ transitions are dominated by the 


distance dynamics (see e.g. 
Assuming that both Kg 


Ref. 20 


7r‘’e'''e 


GP-violating part of the Zsd effective amplitude, we can easily relate their widths. 
Indeed, neglecting the electron mass and the effects of the small vector coupling 
of the electrons to the Z, leads to T{Kl —>• 7r°e'''e“) = T{Kl —> 7r°z/h)/6. Using 
this approximate relation we have derived the upper bound for B{Kl —>■ 7r‘^e’''e“) 
reported in the last line of Table 2. As it is well known, in addition to the direct 
GP-violating transition, the Ki —> 7r‘^e'''e“ decay can proceed through indirect CP- 
violation [Kg Kg —>■ 7T^e'^e~) or via the GP-conserving two-photon exchange 
{Ki —> 7r°77 —>■ 7T^e'^e~). However, both these mechanisms are expected to produce 


corrections to B{Kl 


n^e'^e 


) at the level of fewxlO at most 


This 


ensures that a detection of B{Kl n^e'^e ) above 10 can be considered as a 
clear signature of new-physics. 
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As can be expected, the upper limits for the supersymmetric branching ratios 
shown in Table are much larger than those reported in , where the supersym¬ 
metric contributions to K ^ ttvv have been evaluated essentially without allowing 
(d, s)l — iji mixing. We stress, however, that our upper bounds are signihcantly 
larger also than those recently obtained in [Q, where the (d, s)l — mixing has been 
evaluated only to hrst order within the mass-insertion approximation. 

To conclude this section we emphasize that in these numerical results we did 
not take into account other possible sources of enhancement. Indeed we could have 
obtained even bigger effects playing with the various supersymmetric mass ratios 
(that we have set to 1 just to simplify our results). In particular, larger effects are 
obtained with a wino mass lighter than the average squark mass. Moreover, we have 
neglected possible constructive interferences between the leading (LR)^ terms and the 
subleading, but still not negligible, LR terms. Finally, we have neglected possible 
destructive interferences between chargino- and gluino-mediated amplitudes when 
evaluating the bounds on the LR couplings induced by mixing: this effect 

could easily lead to overcome the stringent constraint in Eq. ( p.l5|) . 


5. Conclusions 

In this paper we have analyzed supersymmetric contributions to rare K decays medi¬ 
ated by an effective Zsd vertex. We have adopted the strategy of the so-called mass- 
insertion approximation, which consists in assuming that the squark mass matrices 
are almost diagonal, and that their diagonalization can be performed perturbatively. 
While recent similar analyses have stopped this approximation to the hrst order, we 
have argued that in the present case it is necessary to go up to the second order in 
this expansion to account for all possible important effects. 

This result does not contradict the validity of the mass-insertion approximation. 
Rather, we have stressed the fact that there is an interplay between the squark mass 
matrices and other mass matrices present in the theory. The reason why the second- 
order terms in this expansion can be more important than the hrst-order ones, is 
because they do not contain anymore off-diagonal CKM matrix elements which are 
known to be suppressed. In other words, we could say that all mass matrices (both 
those of the quarks and of the squarks) in the supersymmetric theory are almost 
diagonal, and that for all these matrices we count off-diagonal elements as of order 
e. According to this counting rule, both the SM and the SUSY contributions to this 
process are of order e^, and here we have for the hrst time presented a complete 
result of SUSY effects at order e^. 

Moreover, for reasons related to the necessary presence of S'f/(2)i-breaking ef¬ 
fects in the effective Zsd vertex, the supersymmetric contributions generated at the 
hrst order in the mass-insertion always appear suppressed by oh-diagonal elements 
of the chargino mass matrix. These vanish as 0{mw/Mg) in the limit of a large 
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supersymmetric scale, Ms, thus providing an additional damping factor which can 
be avoided only at the second order in the expansion of the squark mass matrices. 
This suppression as well as the CKM one can only be avoided considering a double 
LR mixing in the up-squark sector. 

We have performed a numerical analysis of the present bounds on the off- 
diagonal LR elements of the up-squark mass matrix relevant to the effective Zsd 
vertex. As a result, we have found that to our present knowledge the term which 
had been neglected so far (i.e. the one generated at second order in the mass- 
insertion approximation) is the most dangerous one, and could lead to very large 
enhancements in rare kaon decay rates. We have shown that the rate 

could be enhanced up to one order of magnitude with respect to the SM prediction, 
whereas the neutral decay mode ^ could be enhanced by up to two orders 
of magnitude. The same two orders of magnitude enhancement could be produced 
also in the decay Kl ^ Finally, we have also briefly discussed why the 

supersymmetric box contributions to these decays can be neglected as long as we are 
interested in potentially large effects. 

Our results show that the current experimental efforts in the search for these 
rare decays are very much welcome and could give us valuable information on the 
flavor structure of the soft-breaking terms of a generic supersymmetric extension of 
the SM. Interestingly, we will not have to wait too long before experiments will reach 
the sensitivity necessary to observe, or at least to constrain, these supersymmetric 
effects. Indeed a preliminary evidence of the > vr+z/P decay has been recently 

obtained |]^ and the BNL-E787 Collaboration is already analyzing new data on 
this mode. A sensitivity on B{Kl —>■ 7T^e~^e~^) at the level of 10“^° is expected in few 
years by the KTeV experiment at Fermilab [^. Finally, concerning the challenging 
Ki ^ 7r°z/P channel, while waiting for the dedicated experiments aiming to reach 
a sensitivity of 10“^^ even a non-dedicated experiment like KLOE [^] has a 


chance to give new and valuable information on possible extensions of the SM, since 
it can reach a sensitivity of 10“® 
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Appendix 


Expansion of the mass matrices around the diagonal. 

Here we report the formulae needed to make the expansion around the diagonal of 
the mass matrices up to second order, i.e. including two mass insertions. Given an 
n X n Hermitian matrix M, we can decompose it in the form 

M = + M\ (A.l) 


where = diag(m5’,..., m^) and has no elements on the diagonal. M can be 
diagonalized by a unitary matrix X, such that XMX'I' = diag(mi,... ,m„). Then, 
if / is an arbitrary function, we have 


Xlkfijnk)Xkj = Sijf{m°) 

+ , (A.2) 

where we have adopted the notation of Buras, Romanino and Silvestrini to dehne 
an n-argument function from an ? 7 , — 1-argument one: 


f{x,y,Zi,...,Zn-2) 


f{x,zi,...,Zn-2) - f{y, ^ 1 ,■ ■ ■, Zn-2) 

x-y 


(A.3) 


Loop functions. 

The loop functions appearing in the top-quark penguin diagrams discussed in sec¬ 
tion 1^ are given by 


C{x) = 


X 


X 


H{x) = — 


X — 6 3x + 2 ' 

f logo: ^ 1 \ 


{x — ly 


X 


(A.4) 

(A.5) 


The multi-variables functions k{xi,, Xn), j{xi,..., Xn) and Z(xi,..., a;„), occur¬ 
ring in chargino-squark diagrams, are dehned according to the recursive formula 
given in (|A.3|). The explicit expression of the single-variable functions are 


3{x) = 


X logx 
X — 1 


k{x) = X j{x) , 


(A.6) 


. (A.7) 

\X xj X \X xj \X X J \XulI Xu^l J 


18 













References 


[ 1 ] 


[ 2 ] 


[3] 

[4] 

[5] 

[ 6 ] 

[7] 

[ 8 ] 
[9] 

[ 10 ] 

[ 11 ] 


For excellent reviews see e.g.: 

H. Nilles, Phys. Rep. 110 (1984) 1; 

H. Haber and G. Kane, Phys. Rep. 117 (1985) 75; 

R. Barbieri, Riv. Nuovo Cimento 11 (1988) 1; 

S. P. Martin, |hep-ph/9709356| . 


S. Dimopoulos and H. Georgi, Nucl. Phys. B193 (1981) 150; 

J. Ellis and D.V. Nanopoulos, Phys. Lett. BllO (1982) 44; 

R. Barbieri and R. Gatto, Phys. Lett. BllO (1982) 211; 

M.J. Duncan, Nucl. Phys. B221 (1983) 285; 

J.F. Donoghue, H.P. Nilles and D. Wyler, Phys. Lett. B128 (1983) 55. 


L. J. Hall, V.A. Kostelecky and S. Rabi, Nucl. Phys. 267 (1986) 415. 

J.S. Hagelin, S. Kelley and T. Tanaka, Nucl. Phys. B415 (1994) 293; 

E. Gabrielli, A. Masiero and L. Silvestrini, Phys. Lett. B374 (1996) 80; 

F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B477 (1996) 321. 

M. Misiak, S. Pokorski and J. Rosiek, jhep-ph / 9703442] , to appear in Heavy Flavours 
//, ed. A.J. Buras and M. Lindner (World Scientific, Singapore). 


J.A. Bagger, K.T. Matchev and R.-J. Zhang, Phys. Lett. B412 (1997) 77; 
M. Giuchini et al, TUM-HEP-320/98, hep-ph /9808328| . 


Y. Nir and M.P. Worah, Phys. Lett. B423 (1998) 326. 


A.J. Buras, A. Romanino and L. Silvestrini, Nucl. Phys. B520 (1998) 3. 

N. Gabibbo, Phys. Rev. Lett. 10 (1963) 531; 

M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 


T. Inami and C.S. Lim, Prog. Theor. Phys. 65 (1981) 297. 

G. Buchalla and A.J. Buras, Nucl. Phys. B400 (1993) 225; Nucl. Phys. B412 (1994) 
106. 


[12] Y. Nir and D. Silverman, Phys. Rev. D42 (1990) 1477. 

[13] J.A. Gasas and S. Dimopoulos, Phys. Lett. B387 (1996) 107. 

[14] Y. Grossman and Y. Nir, Phys. Lett. B398 (1997) 163. 

[15] G. D’Ambrosio, G. Isidori and J. Portoles, Phys. Lett. B423 (1998) 385. 

[16] S. Adler et al. (E787 Collab.), Phys. Rev. Lett. 79 (1997) 2204. 

[17] V. Kaplunovsky and J. Louis, Phys. Lett. B306 (1993) 269; 

J. Louis and Y. Nir, Nucl. Phys. B447 (1995) 18. 


19 










[18] G. Buchalla and A.J. Buras, Phys. Rev. D54 (1996) 6782. 


[19] A.J. Buras and R. Fleischer, TUM-HEP-275-97, |hep-ph/9704376| , to appear in Heavy 
Flavours II, ed. A.J. Buras and M. Linder (World Scientific, 1997). 


[20] G. D’Ambrosio, G. Ecker, G. Isidori, J. Portoles, hep-ph /9808289| , INFN-NA-IV-98- 
25, to appear in JHEP. 


[21] S. Bertolini and A. Masiero, Phys. Lett. B174 (1986) 343; 
G. Giudice, Z. Phys. C34 (1987) 57; 

1. Bigi and F. Gabbiani, Nucl. Phys. B367 (1991) 3; 

G. Couture and H. Konig, Z. Phys. C69 (1995) 167. 


[ 22 ] 

[23] 


R. Kessler, talk given at the XXVI SLAG Summer Institute, August 3-14, 1998. 


I.-H Chiang et al, BNL-P926 (1996); 
T. Inagaki et al., KEK-96/13 (1996); 
K. Arisaka et al., hep-ex/9709026| . 


[24] The KLOE Collaboration, The KLOE Detector, Technical Proposal, LNF-93/002 
(1993); Status of the KLOE Experiment, LNE-97/033 (1997). 


[25] 


F. Bossi, G. Colangelo and G. Isidori, LNF-98-004-P, |hep-ph/9802345 
EPJC. 


to appear in 


20 













